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Abstract 

CN We perform a chiral extrapolation of the baryon octet and decuplet masses in a relativistic 

O formulation of chiral perturbation theory. A partial summation is assumed as implied by the use of 

physical baryon and meson masses in the one-loop diagrams. Upon a chiral expansion, our results 
Ph are consistent with strict chiral perturbation theory at the next-to-next-to-next-to-leading order. 

All counter terms are correlated by a large-Ac operator analysis. Our results are confronted with 

recent results of unquenched three-flavor lattice simulations. We adjust the parameter set to the 
p ^ pion-mass dependence of the nucleon and omega masses as computed by the BMW Collaboration 

and predict the pion-mass dependence of the remaining baryon octet and decuplet states. The 

^ current lattice simulations can be described accurately and smoothly up to pion masses of about 

> 

QQ 600 MeV. In particular, we recover the recent results of HSC without any further adjustments. 
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I. INTRODUCTION 



QCD lattice simulations offer the unique opportunity to determine low-energy parame- 
ters of the chiral Lagrangian. Currently, various lattice groups work on the baryon octet 
and decuplet masses in unquenched simulations with three light quarks [1-6]. Since the 
simulations are performed also at quark masses larger than those needed to reproduce the 
physical hadron masses, new information is generated that may be used to determine so-far 
unknown low-energy constants. 

The strict chiral expansion of the baryon masses based on the flavor SU(3) heavy-baryon 
formulation is poorly convergent [7-14]. This is in contrast to the two-flavor formulation, 
which appears to be sufficiently well converging as to justify its direct application to lattice 
data [15, 16]. Thus, an extrapolation of the recent lattice simulation [1-5] down to the 
physical limit is not straightforward and requires detailed studies. Indeed, at next-to-next- 
to-leading order (NNLO) the heavy-baryon formulation cannot describe the latest LHP [2] 
and PACS-CS [ '>, 17] lattice data. 

A phenomenological remedy of the three-flavor convergence problem was suggested long 
ago by Donoghue and Holstein [8, 9, 12, 18]. It was shown that the introduction of an 
ultraviolet cutoff of the order of the kaon mass tames the large contribution of the one- 
loop effects. For a recent application of such a scheme to the chiral extrapolation of the 
baryon masses, we refer to Ref. [i ,■]. In our previous work [ ], we suggested a possible 
alternative. A partial summation scheme was suggested where physical meson and baryon 
masses are used in the one-loop expression constituting the NNLO effects. This leads to 
very reasonable results for physical quark masses. In a follow-up work [ ], such a scheme 
was applied to a chiral extrapolation of the baryon octet and decuplet masses. Though the 
overall quark-mass dependence of the early lattice results of MILC [ ' ] was recovered roughly, 
the self-consistency of this approach leads to a striking and unexpected phenomenon. The 
system of eight nonlinear coupled equations implies not necessarily a continuous quark-mass 
dependence of the baryon masses. In the more recent work [21], yet a different strategy was 
investigated. It was shown that with a phenomenological adjustment of the meson-baryon 
coupling constants in the one-loop baryon octet and decuplet self-energies evaluated at the 
NNLO level, one may get close to the simulation results of PACS-CS [ ]. In a strict chiral 
expansion, the effect of modifying the meson-baryon coupling constants off their chiral SU(3) 
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symmetric values enters at next-to-next-to-next-to-next-to-leading order (N^LO). It is yet to 
be clear why the N^LO effects are possibly more important than the next-to-next-to-next- 
to-leading order effects. 

It is the purpose of the present study to extend our previous work [ ] and confront it with 
recent unquenched lattice simulations. We will focus on the results of the BMW group [5] 
since they provide results that can be used directly in the continuum with negligible lattice 
effects. For a recent detailed study of lattice volume effects of the results of the NPLQCD 
Collaboration [22], we refer to Ref. [23]. The incorporation of the next-to- next-to-next-to- 
leading order effects in our approach will allow for a quantitative extrapolation of lattice 
results. At this order, the meson-baryon coupling constants are SU(3) flavor symmetric. 
AUtogether, there are 34 additional counter terms to be considered when going from the 
NNLO to the N^LO level. While the relevant counter terms were constructed previously 
in Ref. [7] for the baryon octet, analogous terms for the baryon decuplet will be presented 
fully for the first time in this work. Some terms relevant in the decuplet sector can be found 
in Refs. [15, 24-27]. Because of the large number of unknown parameters an analysis at 
N^LO appears futile at first. However, we will demonstrate that using large- A^^c correlations 
of the counter terms [26-29], a meaningful extrapolation can be performed. The number of 
additional unknown parameters is reduced significantly. 

At subleading order in the large- A,, expansion, we find the relevance of alltogether 20 
parameters, where we find a subset of 15 parameters to be most important. The latter are 
adjusted to reproduce the empirical baryon octet and decuplet masses together with the 
pion-mass dependence of the nucleon and omega mass as predicted by recent QCD lattice 
simulations of the BMW group [ ]. The smooth pion-mass dependence can be reproduced 
accurately. A prediction for the pion-mass dependence of the remaining octet and decuplet 
masses is presented and confronted with available unquenched three-flavor simulations of 
various lattice groups. 

II. CHIRAL LAGRANGIAN WITH BARYON OCTET AND DECUPLET FIELDS 

The construction rules for the chiral SU{3) Lagrangian density are explained in Refs. 
[7, 13, 30-34]. We recall the terms relevant for our work using the notation and conventions 
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of Refs. [11, 27]. The leading order chiral Lagrangian is of the order Q and is given by 

= tr|5(2^ -M[8])5}+Ftr|fi7^75[zf/^,5]}+Dtr|57'^75{zf/^,5}} 

- tr [b^ -{{lip - M[io]) g^"" - I {YD-" + + Y{^ P + M[,o])Y) 5.} 

+ C (tT{{B^-tU^)B} + R.c}j + HtT^iB^ -^,^5 B^)tU'} , (1) 

with Ufj, = id^^/{2f) + ■ ■ ■ . The leading order baryon masses are given by M[8] and M[io] 
for the members of the flavor SU{3) octet and decuplet, respectively. The parameters F 
and D in Eq. (1) may be determined from the study of semileptonic decays of baryons, 
B ^ B' + e + Ue- This leads to F ~ 0.45 and D ~ 0.80 [35, 36], the values used in this work. 
The value of C may be extracted from the hadronic decays of the decuplet baryons. We 
recall from Ref. [ Ui] the empirical value C = 1.6. The parameter H is poorly determined 
by experimental data so far [37]. Using large- A'^c sum rules, the parameters C and H may 
be also estimated given the empirical values for F and D [ ]. It holds 

H = 9F-3D, C = 2D, (2) 

at sub leading order in the 1/Nc expansion. 

A complete list of chiral symmetry-conserving counter terms, relevant for the calcu- 
lation of the N'^LO baryon mass corrections, was given in Refs. [27, 36]. Following these 
works, we group the counter terms according to their Dirac structure and display the 
relevant terms only. It holds: 

£(^) = £(^) + £(^), (3) 

with 

= tr{5 5}tr|[/^[/'^}-^(7f)tr|s[/'^}tr|f/^5} 

- \ 9^n^tT {b {{U„ U'^},B}]-\ tr {b [{U„ U^} , B] } 
+ ^ /if) tr ■ fi^} tr t/'^} + ^ /if) tr {fi^ ■ 5^} tr jf/'^ f/,} 
+ /if tr I ■ B^) (U" t/.) } + ^ /^f tr I (fi^ ■ B-"^ {[/^ f/,}} 
+ /if)tr|(5^-f/,) (t/'^-i?^)} 

+ ^ /if tr{ ■ U^) (U" ■ 5.) + (fi^ ■ f/^) (U^ ■ fi.) } , 
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- ^f^r^ (tr|5f/^}z7'^tr|f/,D'^5} +tr|5f/,}2 7^tr{f/^D'^5} +H.C. 

- g (^g^^ (tr |5 z Y {{U^, U,} , Z^'^i?} } + H.c. 

- ^ g^P (tr |5z7^ [{f/^, f/.} , D^B] } + H.c.^ 
+ ^-hP (tr|5A-2 7^^'^5^}tr{f/^f/.}+H.c.' 
+ ^ hP (tr I (5, ■ z 7^ D'^i?^) {[/^, t/,}} + H.c. 
+ ^ /if) (tr I (5a ■ f/^) 2 7^ (t/. ■ ^"i?") + [Bx ■ U,)tr [U, ■ D'^B^) } + H.c). (4) 

The large number of unknown chiral parameters at this order is reduced by matching the low- 
energy and the 1/Nc expansions of the product of two axial- vector quark currents [27, 36]. 
The 17 parameters in Eq. (4) are correlated by the 12 sum rules 

hr' = lsr. (5) 

leaving only the five unknown parameters, go^\ gi^\ g^^'' dP idP ■ 

It remains to detail the explicit symmetry-breaking terms. We collect the terms relevant 
for the computation of the baryon self-energies at N^LO. There are five terms: 

= 2 60 tr [B B) tr + 2 6^ tr [B B}) + 2 6^ tr [B B]) 

- 2 rfo tr {B^ ■ B^) tr(x+) - 2 rf^ tr ( {B^ ■ B^) x+) , 
1 
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X+ = Xo - ^ Xo} + O ($^) , (6) 



with Xo = 2 i?o diag(m, m, rrj-s) proportional to the quark-mass matrix. We do not consider 
isospin- violating effects in this work, and we use / = 92.4 MeV. A matching of the chiral 
interaction terms (6) to the large- A'c operator analysis for the baryon masses in Refs. [26, 28] 
leads to one sum rule: 

hD + hF=^dD, (7) 



accurate to subleading order in the 1/Nc expansion. 

We continue with the symmetry-breaking part of the chiral Lagrangian. There are five 
terms and 12 terms: 

= CotT{B (tip - M[s]) B)tT{x+) + CDtT{B {lip - M[s])[B,x+]) 
+ CFtT{B{zl/) -M[s]){B,x+}) 

- ^,tT{B^{t]p -M[,o])B>^)tT{x+)-^Dtr{B^itJp -M[,o])B'^X+), 
Cf = CotT{BB)tT{xl)+c,tT{Bx+)tT{x+B) 

+ C2tT{B{xl,B})+cstT{B[xl,B]) 

+ C4 tr {B B}) tr(x+) + C5 tr {5 [x+, B]) tr(x+) 

+ c,tT{BB) (tr(x+))' 

- eotr {B, ■ B^) tr {xl) - ei tr {{B, ■ x+) (x+ • B^)) 

- 62 tr ( [B, ■ B^) ■ xl) - 63 tr ( [B^ ■ B^) ■ x+) tr 

- e4tr(5^-5'^)(tr(x+))' . (8) 

We consider again large- A^'c sum rules for the parameters introduced in Eq. (8). For Co, (d, Cf, 
and ^o,C,D, the sum rules are analogous to the ones for the terms, i.e. it holds: 

Cd + Cf = 1^d, (9) 

at subleading order. A matching of the chiral interaction terms (8) to the large- A^^c operator 
analysis for the baryon masses in Ref. [2i')] leads to the seven sum rules: 

- 1 - 3 _n 

Co — -Ci, C2 — — - Ci, C3 — u , 

eo = 0, ei = -2c2, 62 = 802, 63 = 8(04 + 05), (10) 

valid at NNLO in the expansion. Assuming the approximate validity of Eq. (10), it suffices 

to determine the five parameters 01,04,05,06, and 04. 
(3) 

The terms in £^ are redundant upon a suitable redefinition of the baryon fields. Why do 
we consider such terms at all? The reason is that the redundance of the parameters Co,Cd,Cf, 
and C,o,^D is lifted once we insist on the large- A'^c relations (10). This is seen by eliminating 
in application of the equation of motion for the baryon fields. A renormalization of 
terms already present in the chiral Lagrangian, in particular, those in Cx \ arises. We find 
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the renormalized or effective coupling strengths: 

cf = -2 (Cd bo - Cf hp) + Co, ef = eo, 

cf = -| (Cd &d - 3 Cf hp) + ci, ef = | + ei, 

cf = 2 (Cd - 3 Cf hp) + C2, ef = -2 ^d + 62, 

cf = -2 (Cd ^f + Cf &d) + C3, ef = -2 (^o c/d + in do) + es, 

cf = -2 (Co 6d + Cd &o + 2 Cd &D - 2 Cf &f) + C4, ef = -2^0 do + e^, 

cf = -2 (Co6^ + C^6o) + C5, 

cf = -2 (Co bo - Cd &D + Cf &f) + ce- (11) 

Inserting Eq. (11) into the sum rules (10) reveals that there is no nontrivial way to dial the 
parameters (dXf, and as to be compatible with Eq. (10). It follows that those param- 
eters are independent of the five parameters ci,C4,C5,C6 and e^. For the singlet parameters 
one finds the correlation 

S Co (bn + bp) =CodD, (12) 

which is consistent with the trivial solution Co = = ^q. The results (11, 12) illustrate that 
the N^LO effect of a variation of the parameters Co and ^0 as correlated by Eq. (12) can be 
reproduced by a suitable variation of the parameters €4^5^ and e4. We conclude that given 
the sum rules, (9, 10) there are alltogether 8 symmetry-breaking parameters at N^LO: the 
parameters Ci, C4, C5, Ce, e4, and Cf,^d together with either ^0 or Co- 
in anticipation of our results, we state the crucial importance of the parameters Cfj^d- 
Within the self-consistent approach applied in this work, they lead to a smooth chiral ex- 
trapolation of the baryon masses. Any attempt with Cf = 'Cd = to establish a smooth 
chiral extrapolation would be futile. 



III. CHIRAL LOOP EXPANSION OF THE BARYON MASSES 

We turn to the computation of the baryon masses. The baryon self-energy, Tib{]^), may 
be considered to be a function of Pfj,'~f^ only, with the 4-momentum of the baryon B. This 
is obvious for the spin one-half baryons, but less immediate for the spin-three-half baryons. 
We refer to Ref. [J J] for technical details. To order Q^, the self-energy receives contributions 
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from tree-level diagrams and one-loop diagrams: 

Sb(Mb) = S*^^^-'^™' + , (13) 

where the index B stands for the members of the flavor SU{3) octet and decuplet, 
Be [8], [10]. The separation of the baryon self-energies into a loop and a tree-level contri- 
bution is not unique depending on the renormalization scheme. In this work, we apply the 
XMS scheme developed in Ref. [I I], where we keep the dependence on the ultraviolet renor- 
malization scale only. A possible dependence on an infrared renormalization scale is ignored 
in this work as to be close to more conventional renormalization schemes. A matching with 
alternative renormalization schemes is most economically performed by a direct comparison 
with the explicit expressions of this section. Our renormalized tree-level self-energies are 
collected in the Appendix. 

The physical mass of the baryon Mb is determined by the condition 

f Mrgi for B e [8] 

[ M[io] for B e [10] 
where M[8] and M[io] are the renormalized and scale-independent bare masses of the baryon 
octet and decuplet. We consistently use a bar for renormalized quantities throughout this 
work. At NNLO, the one-loop contributions probe the coupling constants F,D,C, and H 
introduced in Eq. (1). Complete expressions for baryon octet and decuplet states were first 
established in Ref. [li]. Partial results are documented in Ref. [39]. Here, we complement 
our previous result by additional contributions from the counter terms (4, 8), which turn 
relevant at N^LO. For previous N'^LO studies in the baryon octet sector see Refs. [7, 40]. 
Some partial results in the baryon decuplet sector can be found in [ ]. Alltogether, we 
obtain for the renormalized loop contribution the expressions 

Woop _ V- /cSV/ {Mb + M^, (, , Iq \ , Mj-Ml - ' 

Qe[8],iJe[io] \ / L ^ R Q 

'{Mr - Mb) {Mr + Mb? + m% h Ml + Q MrMr - 2 Ml 



12 Mb Ml 12 Mb Ml 




+ 7^E f^B^--Q^S-I-Q^BGi3)/g, (15) 
^ ■' ' Qe[8] 



and 




+ ' 24MI 24MI ' 



Qe[8],iJe[io] 



+ M|-m|J + 36 M| M| ~ 

{MB + Mnfrn'^Q 3 M| - 2 M| M^? + 3 M| M| - 2 M| 2 
^ 36M|M2 + 36M| Ml 

+ E (^B^ - - ^ ^22] ^ , (16) 
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2M2 / \Ml 

M I-2pqrMb \ _ ( MI + 2pqrMb 
ml + Ml ) rnl + Ml 

4 2 + ^M| ' E^-M^+Pq^. (17) 

The sums in Eqs. (15, 16) extend over the intermediate Goldstone bosons {Q G [8]) baryon 
octet {R G [8]) and decuplet states (i? G [10]). The couphng constants are determined 
by the parameters F, D,C, H . They are hsted in Ref. [11]. The couphng constants 
probe the renormahzed symmetry-breaking parameters bQ,b£,,bF,dQ, and dr,. They are 
detailed in Table I, together with Gq]^ and which are proportional to the symmetry- 
preserving parameters introduced in Eq. (4). In Table I, we apply the notation 

fcf-ftP' + jAf, ftf + hf'^hf> + \hf>. (18) 

The 6 parameters hi_Q enter the decuplet self-energy in the three combinations (18) only. 

The mesonic tadpole Iq has a logarithmic dependence on the ultraviolet renormalization 
scale /iuv- The ultraviolet scale dependence of Eqs. (15, 16) is counteracted by the to-be- 
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TABLE I: Coefficients G^-^' G^^> and G 
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Empirical value at fiuv = 0.8 GeV 


2 Lq — 


-0.1 X 10-3 


2 Lg — L5 


+0.4 X 10-3 


3 L7 + Lg 


-0.3 X 10-3 


m Bq 


9.91 X 10-3 GeV2 


rus Bo 


237xlO-3GeV2 


f 


92.4 MeV 



TABLE II: Low-energy coupling constants used in this work. 



specified tree-level self-energy of Eq. (13). Applying a further chiral expansion to Eq. (13), 
it was demonstrated that the physical masses are renormalization-scale-independent. Here, 
we generalized that result to order N3L0. 

In the loop expressions (15, 16), we use the meson masses accurate to the NNLO as 
derived in Refs. [42, 43]. We recall 
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(19) 



with the renormalized mesonic tadpole integrals It^^Ik-, Ir) as given in Eq. (17). The empirical 
values for the three relevant combinations 2Lq — L^^, 2 Lg — L5, 3 L7 + Lg are given in Table II 
at the renormalization scale ^jjy = 800 MeV. 

We provide a more specific discussion of the renormalized baryon mass parameters that 
enter the mass equation (14). With Mjgj and M[io], we denote the renormalized form of the 
bare parameters M[g] and M[io]. They do not coincide with the chiral SU(3) limit of the 
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baryon masses. The latter are determined by a set of nonlinear and coupled equations: 



where we identified the baryon octet and decuplet masses in the one-loop self-energy with 
M and M + A, respectively. The baryon masses receive contributions from the renormalized 
tree-level terms (14) and from the renormalized one-loop self-energies (15, 16). The sum 
of both provides the scale-invariant chiral limit values of the octet and decuplet masses, M 
and M + A. Given any values for Mjg] and M[io], the chiral mass parameters M and A are 
obtained by a numerical solution of Eq. (20). Though in the chiral limit the self-consistency 
condition (14) is not a significant effect, a perturbative expansion of Eq. (20) in A/M is 
rapidly converging, this changes as we turn on fiavor-breaking effects and increase the quark 
masses. 

We emphasize again that Eqs. (15, 16) depend on the physical meson and baryon masses 
niQ and Mr. This defines a self-consistent summation since the masses of the intermediate 
baryon states in Eqs. (15) and (16) should match the total masses. The baryon masses are 
a solution of a set of eight coupled and nonlinear equations in the present scheme. This is 
a consequence of self-consistency imposed on the partial summation approach. The latter 
is a crucial requirement since the loop functions depend sensitively on the precise values of 
the baryon masses. 

We affirm that a strict chiral expansion of the expressions (15, 16, Al, A2) to N^LO leads 
to results that are renormalization scale-independent. 

IV. QUARK-MASS DEPENDENCE OF THE BARYON MASSES 

We discuss the determination of the parameter set. In this work, we introduced altogether 
41 parameters. At LO there are 2 parameters M[8] and M[io]- At next-to-leading order, there 




(20) 
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are the 5 parameters 6o,d,f and do^D, together with / = 92.4 MeV. The parameters F, D, C, 
and H turn relevant at NNLO. We use the large- A^^^ sum rules (2) together with F = 0.45 and 
D = 0.8. At N^LO, there are 12+5 symmetry-breaking parameters Co_6, eo_4 and Co,d,f,^o,d 
introduced in Eq. (8) and the 17 symmetry-conserving parameters of Eq. (4). As discussed 
in detail in Sec. II, some of the symmetry-breaking parameters are redundant at N^LO. At 
N^LO, there are altogether 36 relevant parameters. Using large- A'^c relations for the N'^LO 
parameters, the number of parameters was reduced significantly down to 20. A parameter 
reduction by about a factor of two was achieved. Still, 20 is a large number and, in this 
work, we will consider a subset of the most important operators only. In this work, we will 
ignore the role of the five symmetry-conserving N^LO parameters. This leaves us with the 15 
parameters, which we adjust to the physical baryon masses and the quark-mass dependence 
for the nucleon and omega mass as predicted by the BMW Collaboration [ ]. With even more 
accurate and complete QCD lattice results it should be possible to determine the full set of 
large- A'c correlated parameters in the near future. The results from the BMW Collaboration 
are shown in Fig. 2 for three different lattice spacings. The approximate independence of 
the lattice spacing we take as a justification to adjust our parameters without any further 
continuum limit extrapolations. 

Given a set of parameters there is no guarantee for a unique solution of Eq. (14) to 
exist. In particular, there may be a discontinuous quark-mass dependence for the baryon 
masses. This is a consequence of the nonlinearities in our approach as introduced by the 
self-consistency condition. Indeed, various discontinuities in the quark-mass dependence 
of the baryon masses were reported in Ref. [20] based on the chiral Lagrangian relevant 
at NNLO. While at NNLO, it is not possible to avoid such a discontinuous quark-mass 
dependence in our approach; we find that at N^LO, there are parameter sets that do lead to 
a smooth quark-mass dependence of the baryon masses. We observe a necessary condition 
for a smooth extrapolation: 

d ^ , 



. .r (21) 



to hold for all octet and decuplet self-energies. Owing to our self-consistency constraint, 
the condition (21) depends on the physical baryon masses and the parameters Co,d,f,^o,z) 
only. In order to analyze the condition (21) in more depths, we consider the three mass 
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FIG. 1: Chiral extrapolation of the nucleon and omega masses. Lattice data are taken from Ref. 
combinations, 

Ai = |MA+iMs-| (M;v-Mh)-| (Ms*-MA-Mn + M=.) , 
A2 = Mq- Mh* - 2 (M=. - Ms.) + Ms. - Ma , 

A3 = Ms* - Ms - Mh. + Mh , (22) 

studied before in Refs. [26, 28]. As shown in Ref. [26], a strict large- Ai'c expansion of the 
baryon masses at NNLO predicts Ai = A2 = A3 = 0. The merit of Ai and A2 lies in 
their independence of all parameters but Co,d,f,^o,d, and M[8]. The mass combination A3 
has only an additional dependence on c/^ — 3 (6^ + 6/5). These properties are a consequence 
of the self-consistency constraint and the large- A^^c sum rules (10), which we use for the 
renormalized coupling constants at the renormalization scale fiuv = ^[s]- 

Using the empirical values for Ai ~ 3.2 MeV and A2 — —6.1 MeV, together with the 
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Fit 1 


Fit 2 


Fit 3 


Fit 4* 


Fit 5 


Fit 6 


M[8] [GeV] 


0.9138 


0.9178 


0.9203 


0.9111 


0.9159 


0.9186 


M[io] [GeV] 


1.0937 


1.0915 


1.0896 


1.0938 


1.0917 


1.0897 


bo [GeV-i] 


-0.9115 


-0.9234 


-0.9355 


-0.9086 


-0.9201 


-0.9321 


Id [GeV-i] 


0.5729 


0.5978 


0.6205 


0.5674 


0.5919 


0.6141 


bp [GeV-i] 


-0.6322 


-0.6542 


-0.6733 


-0.5880 


-0.6100 


-0.6289 


do [GeV-i] 


-0.2246 


-0.2337 


-0.2432 


-0.2300 


-0.2385 


-0.2478 


Jz) [GeV-i] 


-0.3778 


-0.3691 


-0.3586 


-0.3617 


-0.3541 


-0.3446 


CO [GeV-3] 


0.0177 


0.0022 


-0.0135 


0.0176 


0.0020 


-0.0136 


C4 [GeV-3] 


-0.1035 


-0.2269 


-0.3413 


-0.1659 


-0.2861 


-0.3975 


C5 [GeV-3] 


-0.3683 


-0.3173 


-0.2776 


-0.3320 


-0.2824 


-0.2444 


C6 [GeV-3] 


-1.5436 


-1.5455 


-1.5496 


-1.2366 


-1.2417 


-1.2482 


64 [GeV-3] 


-0.2783 


-0.2847 


-0.2890 


-0.2520 


-0.2605 


-0.2662 


Co [GeV-2] 


1.7373 


1.9201 


2.0962 


1.1279 


1.3120 


1.4885 




0.2582 


0.2408 


0.2244 


0.2848 


0.2670 


0.2505 


Cf [GeV-2] 


-0.1949 


-0.2221 


-0.2504 


-0.2221 


-0.2487 


-0.2769 


e'o [GeV-2] 


1.1999 


1.4051 


1.6083 


1.1964 


1.4026 


1.6061 



TABLE III: The parameters are adjusted to reproduce the empirical values of the baryon octet 
and decuplet masses and the lattice results for quark-mass dependence of the nucleon and omega 
as shown in Fig. 1. The different parameter sets follow with Ad^ = —0.2 GeV^^ and Adn = —0.3 
GeV~^ for the first three and last three fits, respectively. The parameter A^o takes the increasing 
values 0.2, 0.4 and 0.6 in both cases. 

\aige-Nc sum rule Co = 3 (Cf + (d) from Eq. (9), we can determine the parameters (f,d in 
terms of the two singlet parameters Co and .^o- In turn, the condition (21) may be considered 
as a function of those parameters Co and ^o only. We find 

^0 > Co '' ^ 1-1070 + 1.1882 lnM[8] for - 1.5 < Co < 1-5 , (23) 

with M[8] measured in units of GeV. The third mass combination with its empirical value 
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FIG. 2: Pion-mass extrapolation of the baryon octet masses. 

A3 ~ —23.9 MeV leads to the condition 

Co = 3.1398 + 2.5604 InM^g] + 0.5800 A^o + Ac^d ^6.2490 + 2.3801 In M[8]) 

- (^M[io] - M[8]) (^0.9017 + 1.4407 lnM[8] - 1.6120 A^o) , 
eo = + A^o , = dn-3bF-3bn, (24) 

where all parameters are assumed in units of GeV. From the large- A^c sum rule (7) we expect 
Ado = 0. As seen from (24) this would lead to unnaturally large values for (q, at least for 
reasonable choices of Mjg] and M[io]- Since the parameter dn enters at NLO it is justified to 
admit a small Ad^ < 0, as it would arise at the next order in the large- iVc expansion. 
In the following, we assume fixed values for Ado and A^o and adjust the remaining 14 
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^^^ BMW ^#^ PACS 
hfn HSC O exp. 

LHPC extrap. 

fff VIII c; 



0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 0.5 

[GeV'] 

FIG. 3: Pion-mass extrapolation of the baryon decuplet masses. 

parameters to the physical baryon masses and the pion-mass dependence of the nucleon and 
omega masses as predicted by the BMW Collaboration. The results are shown in Fig. 1 
for various choices. We find that the pion-mass dependence of the BMW results can be 
reproduced accurately for any given Ad^ and A^o- The size of the fitted parameters are 
collected in Table 111. While the parameters Ado and A^o cannot be determined from the 
BMW results, the request for natural-size parameters favors Fit 4, for which all parameters 
take a reasonable size. The chiral limit values of the baryon octet and decuplet states, 
M and M -|- A, follow from the solution of the set of nonlinear equations (20). Using the 
parameters of Table 111, we find the ranges 

M = 943.9 ± 1.8 MeV, M + A = 1085.8 ± 1.6 MeV . (25) 



1.5 
1.3 
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In Fig. 2 and 3 we confront our results for the baryon octet and decuplet masses with the 
predictions of various lattice groups. The almost invisible bands in the figure are generated 
by the 6 parameter sets as specified in Table 111. We refrain from incorporating any finite 
lattice effects in our present study, so the comparison in Figs. 2 and 3, is in part, of a 
qualitative nature. The spread in the various lattice simulation results may be taken as an 
indication on the size of different finite lattice effects. Most interesting is the comparison of 
our results with the predictions from HSC [ ], for which one may expect the need of only 
minor lattice corrections. We find it encouraging that our approach appears to recover the 
pion-mass dependence of the unfitted baryon masses of HSC [i] reasonably well. 



V. SUMMARY 



We have studied the pion-mass dependence of the baryon octet and decuplet masses based 
on the chiral Lagrangian truncated at N^LO. The large number of parameters was reduced 
significantly in application of large- A^^c sum rules and therewith allowed for a first meaningful 
analysis of recent QCD lattice results. Altogether, we considered 16 parameters, where we 
ignored the small effects of the 5 symmetry-conserving N^LO parameters. In our analysis, we 
used a covariant form of the chiral Lagrangian and the pertinent loop functions relevant at 
N^LO. Owing to a self-consistency condition, which requires the use of physical masses in the 
one-loop functions, a successful reproduction of the recent results of the BMW Collaboration 
on the nucleon and omega mass was achieved. A smooth quark-mass dependence arose upon 
a suitable choice of the symmetry-breaking N^LO parameters. A prediction for the pion- 
mass dependence of the remaining octet and decuplet masses was presented and confronted 
with available unquenched three-flavor simulations of various lattice groups. We recover the 
recent results of the HSC without any further adjustments. 

With additional lattice data, in particular, on the dependence of the baryon masses on 
the strange quark mass, it should be possible to determine the remaining five symmetry- 
conserving parameters and scrutinize in more depths the reliability of the assumed large- iVc 
sum rules. 
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Appendix A: Tree-level baryon self-energy 



We specify the renormalized tree-level self-energies for the baryon octet and decuplet 
states. There are several contributions to J]*''^^"'^'^^' in our scheme. We express the tree-level 
self-energy in terms of the renormalized coupling constants. It holds: 



itree— level 



-itree— level 



-2 Bo 
-4 So 



-itree— level 



-itree— level 



and 



-1 tree— level 



E tree— level 



tree— level 



-4 Bo {bf {2m + rus) + 6^ (m + m,) + bf (m - m,)) 
Co (2 + ml) + C2 (m^ + m^) + C3 (m^ — m^) j 
Co (2 m ms) + CDim + m^) (f (m - m^)^ ^M^v - M[8]^ , 
bf (2 m + m,) + ^ (m + 2 m^)) 
Co (2 + ml) ^ Ci (m - m^)^ ^ (?^^ + 2 m^)^ 



-4So^ 



-2 5o 
-4 So 



-4So^ 
-2 So 
-4 So 



Co (2m + ms) + ^ Cd (m + 2m,) ) (Ma - M[8] 



6f (2m + m,) + 2 6^^m^ 
Co (2 m^ m^) 2 C2 m^) 
Co (2 m + m,) + 2 Cd m) (^M^ - M^s]) , 
bf {2m + ms) + 6^ (m + m,) - bf (m - m,)) 
-4 S^ (co (2 m^ + ml) + C2 (m^ + ml) - C3 (m^ - m^)) 
—2 So ^Co (2 m ms) Co ("^ + "^s) ~ (m-m,))(MH-M[8i), (Al) 



= -4Sn 



-4So^ 



-2 So 
-4 So 
-4So2 
-2 So 
-4 So 



-4So^ 



rf;;^(2m + m,) + rf^^m) 
Co (2 m^ -|- ml) + m/^ 
^0 (2 m + ms) + m) (m^ - M[io] ) , 



df (2 m + m,) + ^ df (2 m + m,)) 

Co (2 m^ m^) ^ ci (m - m,)^ ^ €2 (2 m^ + m^) j 



3 

Ms* 



Ml 



[10] 



Co(2m + m,) + -CD(2m + m,))( 
df (2m + m,) + ^ df (m + 2m, 
Co (2 m^) ^ Ci (m - m,)^ + ^ ^2 ("^^ + 2 m^)^ , 
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-2 Bo (^eo (2 m + m,) + - (m + 2 m,)) (m^* - M[io]) 
= -4 Bo (df (2 m + m,) + m^) 
-4 So (^eo (2 + ml) + es mfj 
-2 So (eo (2 m + m^) + m,) (^Ma - M[io]) , 

6o + ceBo{2m + rus) , 
If + c^Bq {2m + ms), 
do + 64 Bo {2m + mg), 



b'g = bD + C4Bo{2m + ms), 
= dD + esBo(2m + ms) . 
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